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Abstract 

This paper reports very-long-baseline interferometry observations of the radio-loud broad absorption 
line (BAL) quasar J1020+4320 at 1.7, 2.3, 6.7, and 8.4 GHz using the Japanese VLBI network (JVN) and 
European VLBI network (EVN). The radio morphology is compact with a size of ~ 10 pc. The convex 
radio spectrum is stable over the last decade; an observed peak frequency of 3.2 GHz is equivalent to 
9.5 GHz in the rest frame, suggesting an age of the order of ~ 100 years as a radio source, according 
to an observed correlation between the linear size and peak frequency of compact steep spectrum (CSS) 
and giga-hertz peaked spectrum (GPS) radio sources. A low-frequency radio excess suggests relic of past 
jet activity. J1020+4320 may be one of the quasars with recurrent and short-lived jet activity during a 
BAL-outflowing phase. 

Key words: galaxies: active — galaxies: jets — quasars: absorption lines — radio continuum: galaxies 
- galaxies: quasars: individual (SDSS J102027.20+432056.2) 



1. INTRODUCTION 

Broad absorption line (BAL) quasars arc identified 
in rest-frame ultra violet spectra by broad absorption 
troughs displaced blueward from the corresponding emis- 
sion lines, such as Civ and Mgn (Weymann et al., 1991). 
The blueward displacements, sometimes up to ~ 0.2c, are 
attributed to intervening outflows from active galactic nu- 
clei (AGNs) in our line of sight. The intrinsic percentage 
of quasars with BALs is -20% (Hewett & Foltz, 2003; 
Knigge et al., 2008). Two plausible explanations, an orien- 
tation scheme and an evolutionary scheme, were proposed 



and have been long debated. The orientation scheme pro- 
poses that the outflowing BAL wind is preferentially equa- 
torial (Murray et al., 1995; Proga et al., 2000) and can be 
observed as absorption troughs only when the accretion 
disk is viewed nearly edge-on. The evolutionary scheme 
proposes that BAL outflows are associated with a rela- 
tively short-lived (possibly episodic) evolutionary phase 
(Gregg et al., 2000, 2006). 

Radio observations offer several inclination indicators 
and age estimators that can be used to test these two con- 
troversial schemes for BAL quasars (Stockc et al., 1992; 
Richards et al., 2001; Menou et al., 2001; Brotherton et al., 
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2006; DiPompeo et al., 2012, and references as follows). 
Several radio sources in BAL quasars exhibit rapid ra- 
dio variability that indicates unusually high brightness 
temperatures, which requires Doppler beaming of nearly 
pole-on viewed jets, i.e., accretion disks with small in- 
clinations (Zhou et al., 2006; Ghosh & Punsly, 2007). 
A wide range of spectral indices, including flat spec- 
tra and steep spectra, is consistent with the wide range 
of orientations (Becker et al., 2000; Montenegro-Montes 
et al., 2008; Fine et al., 2011). Although there is weak 
evidence that the spectral indices of BAL quasars are 
steeper than those of non-BAL quasars, which mildly fa- 
vors edge-on orientations (DiPompeo et al., 2011; Bruni 
et al., 2012), a single edge-on geometry cannot describe 
all BAL quasars. Becker et al. (2000) suggested a picture 
in which BAL quasars represent an early stage in the de- 
velopment of quasars on the basis of their compact radio 
morphology, observed in most cases. Montenegro-Montes 
et al. (2008) demonstrated that many radio-emitting BAL 
quasars share several radio properties common to young 
radio sources, such as gigahertz-peaked spectrum (GPS) 
and compact steep spectrum (CSS) sources. GPS sources 
are compact (J$ 1 kpc) and show a convex radio spectrum 
that peaks between 500 MHz and 10 GHz; CSS sources 
are larger (~ 1-20 kpc) and have convex spectra that tend 
to peak at lower frequencies, typically < 500 MHz (O'Dea, 
1998). On the basis of the rarity of Fanaroff- Riley Class II 
radio galaxies in BAL quasars and their observed anticor- 
rclation between the balnicity index and radio loudness, 
Gregg et al. (2006) suggested a model in which a BAL 
phase evolves into a radio-loud phase with a relatively 
short overlap. In this context, AGN wind-induced feed- 
back in the early stages of radio source evolution is dis- 
cussed for the galaxy-black hole coevolution (e.g., Lfpari 
& Terlevich, 2006; Holt et al., 2008). On the other hand, 
Bruni et al. (2012) reported that the fractions of GPS 
candidates are similar in their BAL and non-BAL quasar 
samples, suggesting that BAL quasars are generally not 
younger than non-BAL quasars. Thus, neither the ori- 
entation scheme nor the evolutionary scheme has been 
conclusively demonstrated to date on the basis of these 
inclination indicators and age estimators. 

Very-long-baseline interferometry (VLBI) with milliarc- 
second (mas) angular resolution provides exclusive and 
crucial opportunities for the investigation of the par- 
sec (pc)-scale regions in which the phenomena that these 
inclination indicators and age estimators rely on are ac- 
tually occuring. Several VLBI observations have been 
reported for BAL-quasar radio sources (Jiang & Wang, 
2003; Kunert-Bajraszewska & Marecki, 2007; Kunert- 
Bajraszewska et al., 2010; Liu et al., 2008; Doi et al., 2009; 
Reynolds et al., 2009; Montenegro-Montes et al., 2009; 
Gawronski & Kunert-Bajraszewska, 2010; Yang et al., 
2012; Hayashi et al., submitted) and revealed various sig- 
natures, such as blazar-like jets with a one-sided mor- 
phology, polarized radio emissions, a two-sided morphol- 
ogy suggesting inclined jets, CSS-like characteristics, and 
interactions between interstellar medium and jet. Doi 
et al. (2009) reported their systematic VLBI detection 



survey at 8.4 GHz for 22 radio-loud BAL quasars us- 
ing the Optically ConnecTed Array for VLBI Exploration 
project (OCTAVE: Kawaguchi 2008). The samples of 
bright radio sources (> 100 mJy beam -1 ) were selected 
by position matching between the BAL quasar catalog of 
Trump et al. (2006) from the Sloan Digital Sky Survey 
Third Data Release (SDSS DR3; Abazajian et al., 2005) 
and Faint Images of the Radio Sky at Twenty-cm (FIRST; 
Becker et al., 1995). Most sources (20/22) were detected 
with the OCTAVE baselines, suggesting brightness tem- 
peratures of greater than 10 5 K, and the simultaneous 
coexistence of BAL outflows and nonthermal jets. Four 
sources exhibited inverted spectra, suggesting blazars with 
pole-on-viewed relativistic jets or GPS sources as young 
radio sources. 

In the present paper, we report multifrequency VLBI 
observations of the BAL quasar J1020+4320 (SDSS 
J102027.20+432056.2 at z = 1.962), which was one of 
the four radio sources showing inverted spectra in the 
OCTAVE study. J 1020+4320 exhibits BAL troughs of 
an average velocity of 21341 km s _1 with a width of 
1168 km s ; such a broad width is attributed to an in- 
trinsic absorption. The absorption index (Al; Hall et al., 
2002) is 716 km s _1 , according to the modified defini- 
tion of the Al by Trump et al. (2006). On the other 
hand, the balnicity index (BI; Weymann et al., 1991) 
is (Gibson et al., 2009), according to the criterion of 
a velocity width > 2000 km s _1 in the BI definition. 
J1020+4320 is a compact radio source in the FIRST im- 
age with a resolution of ~ 5". Previous radio observa- 
tions in arcsec/arcmin resolutions (Marecki et al., 1999; 
Vollmer et al., 2008; Orienti et al., 2010; Stanghellini et al., 
2009) suggested that J1020+4320 is a candidate for a high 
frequency pcaker (HFP; Dallacasa ct al., 2000, with a 
spectral peak occurring at frequencies above a few GHz), 
which is considered to be a younger subclass in the GPS 
population. The present paper is organized as follows. 
In Section 2, our observations and data reduction proce- 
dures are described. The observational results are pre- 
sented in Section 3, and their implications are discussed 
in Section 4. Throughout this paper, a ACDM cosmology 
with H = 70.5 km s" 1 Mpc" 1 , fl M = 0.27, and fl A = 0.73 
is adopted (Komatsu et al., 2009). The comoving distance 
is 5225 Mpc; 1 milliarcsecond (mas) corresponds to 8.6 pc 
at the distance from J1020+4320. 

2. OBSERVATIONS AND DATA REDUCTION 

VLBI observations at 1.7, 2.3, 6.7, and 8.4 GHz were 
conducted during a half-year period in 2008 (Table 1). 
J1020+4320 was observed at 1.7 GHz using the European 
VLBI Network (EVN) and the Multi-Element Radio 
Linked Interferometer Network (MERLIN) simultane- 
ously in the snapshot mode. The left- and right-circular- 
polarization signals with a total bandwidth of 32 MHz 
each for the EVN and 16 MHz each for the MERLIN 
were obtained; only Stokes-I correlations were used in 
this study. Observations at 2.3, 6.7 and 8.4 GHz us- 
ing the Japanese VLBI Network (JVN: Fujisawa 2008) 
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Table 1. List of Observations 



V 


Date 


Array 


Antenna 


(GHz) 








(1) 


(2) 


(3) 


(4) 


1.666 


2008Mar02 


MERLIN 


MK LO CA DE KN DA TA 


1.658 


2008Mar02 


EVN 


LO Wb Ef On Mc Tr CA Nt 


2.272 


2008Augl7 


JVN 


VERA x 4 Ud 


6.672 


2008May04 


JVN 


VERA x 4 Ud YM 


8.408 


2008Augl7 


JVN 


VERA x 4 Ud Uc Ks 



Col. (1) center frequency: Col. (2) observation date; Col. (3) array; Col. (4) participating antenna. MK: Mark 2 32 X 25 m, LO: Lovcll 
76 m, CA: Cambridge 32 m, DE: Dcfford 25 m, KN: Knocking 25 m, DA: Darnhall 25 m, TA: Tablcy 25 m, Wb: Westerbork, Ef: 
Effelsberg 100 m, On: Onsala 85 m, Mc: Medicina 32 m, Tr: Torun 32 m, Nt: Noto 32 m, VERA: VLBI Exploration of Radio 
Astrometry (VERA) 20 m, Ud: JAXA Usuda 64 m, YM: Yamaguchi 32 m, Uc: JAXA Uchinoura 34 m, Ks: Kashima 34 m. 



were conducted in single circular polarization at a data 
recording rate of 64 Mbps with 2-bit quantization; this 
provided an observing band-width of 16 MHz for each 
band. The data at 2.3 and 8.4 GHz were obtained simul- 
taneously. Correlation processing was performed using 
the Mitaka FX correlator (Shibata et al., 1998) at the 
National Astronomical Observatory of Japan (NAOJ). 

Data reduction was performed using the Astronomical 
Image Processing System (AIPS; Greisen 2003). For the 
EVN data at 1.7 GHz, amplitude calibration, bandpass 
calibration, flagging, and fringe-fitting were performed in 
the standard manner. The assumed uncertainty of the 
amplitude calibration is 10%. For the JVN data at 2.3, 
6.7, and 8.4 GHz, a-priori amplitude calibration was not 
used because of the lack of a monitoring system of system- 
noise temperature of several antennas at that time. The 
amplitude-gain parameters relative to each antenna were 
obtained from the self-calibration solution for a point-like 
strong source, J0958+4725, which was near the target in 
the sky and was scanned every several tens of minutes to 
monitor the time variation in the system equivalent flux 
density (SEFD) of each antenna. The amplitude scal- 
ing factor was based on the total flux density of OJ 287, 
which was measured with the Very Large Array (VLA) at 
8.4 GHz within a week of the JVN observations and Usuda 
64-m single-dish observations at 2.3, 6.7, and 8.4 GHz 
during the JVN observations. Flux scaling was performed 
by comparing the self-calibration solutions on the scans 
of OJ 287 1 and J0958+4725 at nearly the same eleva- 
tion; a structure model obtained using the Very Long 
Baseline Array (VLBA) was used for the self-calibration 
of OJ 287. Although the structure of OJ 287 is known to 
be variable, it can be assumed to be stable in the range 
of JVN baselines (less than 50 MA at 8.4 GHz). As a 
check, we found that an amplitude scaling factor obtained 
from auto-correlation data for methanol masers, whose 
flux density was determined by a single-dish observation 
with the Yamaguchi 32 m, showed only a 4% difference 

1 OJ 287 is one of the rare objects for which almost all of the 
total flux density can be retrieved at the shortest VLBI base- 
lines according to comparisons between the VLBI correlated flux 
densities and the total flux densities in the VLA and single-dish 
monitoring (UMRAO) for many years. 



from the OJ 287-based scaling factor (Sugiyama et al., 
2011). In addition, the flux scaling factor was approxi- 
mately equal to the factors (with a scattering of less than 
10%) obtained by several other JVN observations (Doi 
et al., 2006, 2007; Tsuboi et al., 2008; Sudou & Edwards, 
2009; Niinuma et al., 2012). We assumed that the uncer- 
tainty of the flux densities for our JVN observations is also 
10%. Imaging was performed using the Difmap software 
(Shepherd et al., 1994) to apply iteratively deconvolution 
and self-calibration procedures. 

3. RESULTS 

3.1. Radio morphology 

J1020+4320 is nearly unresolved in all our VLBI im- 
ages at 1.7, 2.3, 6.7, and 8.4 GHz; the measured to- 
tal flux densities are listed in Table 2. At 8.4 GHz 
(Figure 1), the deconvolved size determined using the 
AIPS task JMFIT in the image domain 2 is 0.6 ±0.2 mas 
at a position angle (PA) of PA = 87° ± 55° and an axis 
ratio of < 0.1, suggesting an elongated structure of ~ 5 pc. 
We also analyzed archival data from the VLBA Imaging 
and Polarimetry Survey (VIPS; Helmboldt et al., 2007) at 
5 GHz and measured the deconvolved size to be 1.1 mas 
elongated at PA =67°, which is consistent with the JVN 
result. The difference in flux density between (simulta- 
neous) MERLIN and EVN observations at 1.7 GHz was 
29 ± 19 mJy, which was minor (~ 15%) as compared to 
the total emission (193 mJy). Thus, most of the emission 
is concentrated in a central region within ~ 25 mas (the 
beam size of EVN) corresponding to ~ 200 pc. An un- 
resolved structure in the MERLIN image constrains the 
entire size to < 190 mas, corresponding to < 1.6 kpc. 

3.2. Radio spectrum 

Our VLBI results suggest a spectrum with a turnover 
(Figure 2). We applied spectral fitting to only our 2008 
VLBI data at four frequencies (see Tables 1 and 2). A 
simple power-law spectral model including synchrotron 
self-absorption (SSA), S v = 5 ^ 2 - 5 [l - exp (-r*/*°- 2 - 5 )], 
provided a peak flux density S p = 417 mJy at a peak 

2 We obtained nearly the same result using the visibility-based 
modelf it in difmap. 
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Table 2. Results of Observations 



V 


Sy 


a 




P.A. 


(GHz) 


(mJy) 


(mJy beam -1 ) 


(mas x mas) 


(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


1.6661 


193 ±19 


0.4 


136.6 x 188.7 


48 


1.6581 


164 ±16 


1.1 


26.6 x 21.3 


-87 


2.272R 


330 ± 34 


3.7 


16.1 x 9.2 


-61 


6.672L 


252 ±26 


2.7 


6.0 x 3.1 


-69 


8.408R 


206 ±21 


1.2 


5.4 x 3.3 


-52 



Col. (1) center frequency. "L" denotes left circular polarization, "R" denotes right circular polarization, and "I" denotes dual circular 
polarization; Col. (2) flux density; Col. (3) image rms noise; Col. (4) Half-power beam width; Col. (5) position angle of major axis of 
beam width. 



frequency v v = 3.2 GHz and ao = — 1.1. A spec- 
tral model including free-free absorption (FFA), S v = 
SqVq exp (— r^ -21 ), provided S p = 392 mJy at v v = 
3.2 GHz 3 and ao = —1.2. Although these spectral fittings 
were unable to discriminate conclusively between SSA and 
FFA, peak frequencies of v v = 3.2 GHz arc suggested in 
either case. 

3.3. Variability 

Our VLBI measurements and the two fitted spectra are 
in good agreement with previous studies (small symbols in 
Figure 2), even those with different beam sizes. Therefore, 
J1020+4320 is a compact and stable radio source without 
significant variability over the last decade. 

4. DISCUSSION 

The observational results confirmed a series of charac- 
teristics corresponding to a young radio source: (1) com- 
pact morphology, (2) a giga-hcrtz peaked spectrum, 
and (3) little variability (O'Dca, 1998, for a review). 
Therefore, we conclude that the BAL quasar J1020±4320 
possesses a young radio source at its nucleus. 

The determined spectral peak frequency of 3.2 GHz in 
the observer frame is equivalent to a peak at 9.5 GHz 
in the rest frame at z — 1.962. The age estimator based 
on the observed correlation between the linear size and 
peak frequency of CSS and GPS radio sources (O'Dea & 
Baum, 1997; Snellen et al., 2000) indicates that the radio 
source in J1020±4320 might be extremely young with an 
age of the order of ~ 100 years (~ 10 pc in size at a ten- 
tative expansion velocity of ~ 0.3c; see O'Dea & Baum 
1997). Indeed, the marginally deconvolved structure with 
an elongation of ~ 1 mas in the JVN 8. 4- GHz and VIPS 
5-GHz images (Section 3) is suggestive of a mini radio 
galaxy of ~ 10 pc across. Alternatively, the result of SSA 
spectral fitting suggests a source diameter of the order 
of ~ 10 pc (with a magnetic field of ~ 0.1 G) under the 



The spectral peak frequencies of J1020+4320 have been deter- 
mined to be ~ 1-5 GHz in a series of previous studies as well 
(Marecki ct al., 1999; Vollmcr et al., 2008; Stanghcllini et al., 
2009; Orienti et al., 2010), which are consistent with our fitting 
results. These minor discrepancies are attributed to their dif- 
ferent model functions, (such as a broken power-law, hyperbola, 
and parabola) which do not have any direct physical significance. 



J1 020+4320 JVN 8.4 GHz 




20 10 -10 -20 

Right Ascension (mas) 

Fig. 1. JVN image at 8.4 GHz. Contour levels arc —1, 
1, 2, 4, 8, 16, and 32 times 3cr of the image rms noise 
(<r = 1.2 mJy beam -1 ). 

near equipartition condition between the energy densities 
of radiating electrons and magnetic field in a homoge- 
neous, self-absorbed, incoherent synchrotron radio source 
with a power-law electron energy distribution. This lin- 
ear size is also consistent with our VLBI images. Thus, we 
can naturally understand the observed radio properties of 
J1020+4320 in terms of its compactness (with a size of 
~ 10 pc). To confirm that the radio source is actually 
young, further investigations might be essential for the 
determination of its kinematic age by, e.g., measuring ex- 
pected spectral changes on spectral property undergoing 
adiabatic expansion (Orienti et al., 2010). 

On the other hand, the data point at 408 MHz (51 ± 
30 mJy with MERLIN; Marecki et al. 1999) is signifi- 
cantly distant from our fitted model spectra. One more 
component is necessary to fit the low-frequency emis- 
sion. It may originate in an optically thin component 
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Fig. 2. Radio spectrum of J1020+4320. Large filled squares 
and a large open square: our VLBI and MERLIN measure- 
ments, respectively, in 2008. Small filled squares: OCTAVE 
at 8.4 GHz on November 04, 2007 (Doi et al., 2009) and 
VIPS at 5 GHz on May 01, 2006 (Helmboldt et al., 2007). 
Small open circles: quasi-simultaneous VLA on June 25, 1999 
(Stanghcllini et al., 2009). Small open triangles: quasi-simul- 
taneous VLA in January, 2004 (Oricnti et al., 2010). Small 
open squares: NVSS on November 15, 1993 (Condon et al., 
1998) and FIRST on February 20, 1997 (Becker et al., 1995) 
at 1.4 GHz and MERLIN at 408 MHz from November 1994 to 
January 1995 (Marecki et al., 1999). Solid and dashed curves 
indicate spectra fitted to our four-frequency VLBI data (large 
filled squares) with spectral models of synchrotron self-ab- 
sorption and free-free absorption, respectively (Section 3). 

at large scales with a steep spectrum, which can also be 
inferred from the observed 29 ±19 mJy at 1.7 GHz (the 
MERLIN-EVN differential flux density) from the region 
of ~ 0.2-1.6 kpc (Section 3). The low- frequency com- 
ponent could be conceivably be explained as a relic of 
past jet activity (Baum et al., 1990; Saikia & Jamrozy, 
2009, and references therein). In the case of J1020+4320, 
the age of the GPS component is ~ 100 years, while the 
low-frequency component may be relic of jets emanated 
~ 0.7-5 x 10 4 years ago (assuming an expansion rate of 
0.1c). The coexistence of components with distinct ages 
indicates recurrent jet activity. The current stage of 3C 84 
may be the closest example for the restarted jet activity 
in the recent 100 years in an existing radio galaxy (Asada 
et al., 2006; Nagai et al, 2009, 2010); the classical double- 
lobcd radio galaxy 1245+676 also contains a GPS compo- 
nent of an expanding compact double with a separation of 
14 pc, suggesting a kinematic age of ~ 200 years (Marecki 
et al., 2003). Similarly, J1020+4320 may not be so young 
as a radio source but contains a very young radio compo- 
nent originating in the restarted jet activity. 

Composite spectra with GPS and MHz components 
have been found in a proportion of radio-loud BAL quasar 
(Montenegro-Montes et al., 2008; Bruni et al., 2012) as 
well as J1020+4320. The fact that the restarted jets 
are associated with BAL features suggests relevance to 
the evolutionary scheme for the origin of BAL quasars 
(Hayashi et al., submitted). In the evolutionary scheme, 
the BAL fraction in the quasar population (^20%; e.g., 



Hcwett & Foltz, 2003; Kniggc et al., 2008) requests the 
duration of the BAL phase of ~ 2 x 10 7 years in the 
AGN lifetime of ~ 10 8 years. On the other hand, time 
scales of ~ 10 6 years for the intermission of the jet have 
been inferred in several double-double radio galaxies such 
as PKS B1545-321 (Safouris et al., 2008), 4C 02.27 
(Jamrozy et al., 2009), B1834+620 (Schoenmakers et al., 
2000), and Cygnus A (Steenbrugge et al, 2008), although 
various different time scales have also been inferred from 
existing observations (Saikia & Jamrozy, 2009, for a re- 
view). Thus, because the time scale between successive 
episodes of the jet activity are possibly much shorter than 
the BAL phase, short-lived GPS sources may appear re- 
peatedly during the BAL phase. J1020+4320 may be one 
of the quasars with a recently reactivated jet during a 
BAL phase. At this stage, our study on a quasar in the 
BAL-GPS composite phase constitutes the first step to- 
ward additional extensive research rather than the con- 
clusive evidence for the relationship between the processes 
initiating/interrupting BAL outflow and a nonthermal jet 
for the galaxy-black hole coevolution. 
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